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ABSTRACT 

We report on a large stellar flare from the nearby dMe flare star EV Lac 
observed by the Swift and Konus-Wind satellites and the Liverpool Telescope. 
It is the first large stellar flare from a dMe flare star to result in a Swift trigger 
based on its hard X-ray intensity. Its peak ix from 0.3-100 keV of 5.3x10"^ 
erg cm^^ s^^ is nearly 7000 times larger than the star's quiescent coronal flux, 
and the change in magnitude in the white filter is >4.7. This flare also caused 
a transient increase in EV Lac's bolometric luminosity (Lboi) during the early 
stages of the flare, with a peak estimated Lx/Lboi ~3.1. We apply flare loop 
hydrodynamic modeling to the plasma parameter temporal changes to derive a 
loop semi-length of l/R-i, = 0.37±0.07. The soft X-ray spectrum of the flare 
reveals evidence of iron Ka emission at 6.4 keV. We model the Ka emission as 
fluorescence from the hot flare source irradiating the photospheric iron, and derive 
loop heights of h/R^=0.1, consistent within factors of a few with the heights 
inferred from hydrodynamic modeling. The Ka emission feature shows variability 
on time scales of ~200 s which is difficult to interpret using the pure fluorescence 
hypothesis. We examine Ka emission produced by collisional ionization from 
accelerated particles, and find parameter values for the spectrum of accelerated 
particles which can accommodate the increased amount of Ka fiux and the lack 
of observed nonthermal emission in the 20-50 keV spectral region. 

Subject headings: stars: activity — stars: coronae — stars: flare — stars: late- 
type — stars: individual (EV Lac) — X-rays: stars - Swift 



1. Introduction 

The origin of solar and stellar flares is thought to be the impulsive release of energy 
from a non-potential magnetic fleld due to an instability arising from magnetic loop foot- 
point motions in the photosphere, resulting in plasma heating, particle acceleration, shocks 
and mass motions, involving the entire outer atmosphere, and producing emissions across the 
electromagnetic spectrum. Flares of all sizes are useful probes of coronal structures, under the 
assumption that there is no difference in the flaring and non-flaring coronal structures during 
the impulsive and transitory energy input. Extreme flare events are useful for revealing the 
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extent to which solar flare physics can be extrapolated (and thus the applicability of the 
solar model in a variety of stellar environments), and for enabling stellar flare diagnostic 
techniques which require high signal-to-noise. 

It has long been known that the largest observed stellar flares on "normal" stars (i.e. 
non-degenerate, non-accreting stars) represent energetic releases which are orders of magni- 
tude more than typical solar flares — up to a factor o f 10^ more energ etic (compared to a 
typic al solar flare energy of 10^^ erg; lEmslie et al.l 120041 ). i.e. ~ 10^* erg (JKiirster &: Schmitt 



19961 ) . Yet the energy radiated as a result of plasma heating usually repr esents only a small 
fraction of the total energy involved in a solar flare: lEmslie et al.l ( 120041 ) estimate that the 
energy contribution of nonthermal particles in solar flares is ~ a factor of 10 hi gher than the 



amount of th ermal plasma heating as expressed in their soft X-ray emission, and lWoods et al. 



(I2OO4J . I2OO6I ) showed that the total solar flare energy can be up to 110 times that seen in the 
GOES soft X-ray band from 1.2-12 keV. 

White light solar flares are difficult to detect in total solar irradiance measurements due 
to the contrast between the flaring emission and total solar disk emission, whereas the red 
nature of the underlying M dwarf spectrum makes the blue flare continuum emission easier 
to detect. In solar white light flares there is a timing and spatial correlation between the 
white light flare emiss ion and signatures of nonthermal emission such as hard X-ray emission 
(JNeidig fc Kanelll993l ). and radiative hydrodynamical simulations of the impact of electron 
beams in the lower atmosphere can reproduce the observed solar continuum enhancements 
(jAUred et al.ll2005[ ). but extension of the same models to flares o n M dwarfs has not to date 
been able to reproduce the observed continuum enhancements (jAUred et al.l 120061 ). Large 
flares from M dwarf stars have been known at optical wavelengths for several decades; the 
UV/optical/infrared spectr al energy distribution ty pically peaks at short wavelengths, with 
a roughly blackbody shape (JHawley fc Fisheiill992l ) having a blackbody temperature Tbb of 
~10'^K. Enh ancements of great er than a few magnitudes are large amplitude events which 
occur rarely; IWelsh et al.l (120071 ) found 7 out of 52 flares in the GALEX data archives to have 
Amuv > 4.4. The most extreme stellar flare ever observed at ground-bas ed near-UV/optical 
wavelengths may be that of V374 Peg (JBatyrshinova &: Ibragimovll200ll ). at AU=11.0. 



Soft X-ray flares result from the transient plasma heating that occurs as a result of energy 
deposition lower in the atmosphere. In the corona, magnetic fleld pressure exceeds the gas 
pressure and therefore dominates the dynamics. The non-flaring coronal emission is easy 
to understand in terms of an optically thin, multi-temperature, and possibly multi-density, 
plasma in coronal equilibrium. The X-ray spectral region also contains diagnostics which 
have been used recently in large stellar f lares to diagnose th e presence of nonthermal electrons 
operating during a large stellar flare (lOsten et al.l 120071 ) as well as probe the interaction 
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between the corona and the lower atmosphere (jOsten et al.ll2007l : lTesta et aLll2008[ ). Changes 
in the characteristics of the flaring plasma with time lead to deductions about heating time- 
scales, place constraints on spatial scales, and provide evidence for temporary enrichment 
of coronal material with plasma from the lower layers of the atmosphere. "Superflares" , 
in which the soft X-ray luminosity increases by > a factor of 100 compared to its normal 



value, are rare but have been seen ori several different active stars (iFavata fc SchmittI Il999 



Franciosini et all I2OOII : Maggio et all l2000t iFavata et all I2OOOI : lOsten et all 120071 ). As the 
quiescent X-ray luminosity of active stars is already at a high level relative to the underlying 
bolometric luminosity (Lx/Lboi ~ 10^^), such an increase during a flare produces a temporary 
but significant increase in the total luminosity of the star. 

The nearby (d=5.06 pc) flare sta r EV Lac has a recor d of frequent and large outbursts 
across the electromagnetic spectrum. iKodaira et al.l ( 119761) reported a long-lasti n g (At > 3 
hrs) flare with a peak enhancement of AU= 5.9, while iRoizman fc Shevchenkol (119821) ob 



served a flare from EV Lac with AU= 6.4 at peak which lasted ~ 9-11 hours. 



Osten et al. 



( I2OO5I ) noted frequ ent soft X-ray flare s from EV Lac. Its bolometric luminosity Lboi has 
been estimated by IFavata et al.l ( I2OOOI ) as 5.25x10^^ erg s~^. Large-scale magnetic flelds 
covering a substantial fracti on of the surface (5 ~ 4 kG, / ~50%) have been measured 
( iJohns-KruU fc Valentilll996l ). The extreme magnetic activity as revealed by the flare obser- 
vations is mirrored in complementary measures of changes in the large-scale magnetic fleld 
configuration on relatively short time-scales of less than an hour (jPhan-Bao et al.ll2006l ). 



Here we describe the properties of a flare serendipitously observed on EV Lac by the 
Swift Gamma-ray Burst satellite, the Konus GRB experiment on the Wind satellite, and 
observed in a fast follow-up by the Liverpool Telescope. The flare is remarkable in the 
magnitude of the peak enhancement at optical and soft and hard X-ray energies. The current 
paper provides an overview of the datasets involved and their reduction (§2), determinations 
of variations in plasma parameters during the flare, accompanying analysis of X-ray spectral 
features, and derivation of parameters from optical data (§3). Section 4 discusses flare energy 
estimates, compares the current flare with other large-scale flares, and the use of the Fe Ka 
6.4 keV line to diagnose conditions in the flaring plasma, while Section 5 summarizes our 
conclusions. 



Data Sets 



An overview of the telescopes and observatories which obtained data for this EV Lac 
flare is given in Figure [H which delineates the timing and relative intensity of the several 
X-ray and 7-ray observations made by Swift and Konus, and also the timing of the UV and 
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optical observations made by Swift and on the ground from the Liverpool Telescope. All 
timing is referenced to the Burst Alert Telescope (BAT) trigger time, 2008-04-25 05:13:57 
UT; we refer to this as TO. 



2.1. Konus-Wind 



Konus (lAptekar et al.lll995l ) is an instrument on the Wind spacecraft which provides 
omni-directional and continuous coverage of cosmic gamma-ray transients. The flare was 
detected in waiting mode by the S2 detector which observes the north ecliptic hemisphere. 
The background is very stable in the S2 detector. There is a statistically significant bump 
in count rate in the lowest energy range (18-70 keV) from ~T0-100:T0+140 s of 5875±692 
counts, which likely indicates the beginning of the superflare. There is also a possible tail 
until ~T0+1000 s. In addition, there may be a precursor event from ~T0-700:T0-400 s 
with 2117±791 counts. Timing analysis shows that the propagation time from the Wind 
spacecraft to the Swift spacecraft is 2.125 s, so the BAT trigger time (UT 05:13:57) corrected 
for time-of-flight is UT 05:13:54.875 at Konus. The peak in the Konus light curve shown in 
Figure [1] indicates that the flare peaked ~60 seconds before the BAT trigger. The Konus 
data can be used to estimate the onset and rise time of the hard X-ray emission from the 
flare, and indicate a relatively fast rise time of less than 2 minutes. 



2.2. Swift/BAT Data 



Swift's Burst Alert Telescope ( iBarthelmy et al.ll2005[ ) triggered on the flare from EV Lac 



during a pre-planned spacecraft slew; the trigger occurred at 2008-04-25 UT 05:13:57 =T0. 
There was a short pointing during which the discovery image was made with BAT, a second 
slew to place the star in the apertures of the narrow-field instruments (X-Ray Telescope 
[XRT] and UV Optical Telescope [UVOT]), and a second longer pointing which overlaps the 
XRT data. The hard X-ray source was already declining in flux as it entered the BAT field of 
view during the pre-planned slew, with a peak intensity in the BAT 15-50 keV band of about 
0.15 photons cm^^ s^^. Comparison of the BAT and Konus/Wind light curves in Figured] 
shows that the BAT triggered on the flare near its maximum in the hard X-ray band. Light 
curves were extracted in the most sensitive 14-30.1 keV range where the spectral distribution 
peaks, and were cleaned of the contaminating sources Cyg-Xl, Cyg-X2, and 3A 0114+650 
which were in the partially coded field of view. Spectra from 14-100 keV were extracted for 
four time intervals, indicated in Figure [H for which there was enough signal. 
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2.3. Swift/XRT Data 



The X-Ray Telescope fJBurrows et al.ll2005l ) started to observe the source at 05:16:48 UT 
(166s after the Swift-BAT trigger, TO). The XRT observes from 0.3-10 keV using CCD 
detectors, with energy resolution in the iron K region (6 keV) of ~ 140 eV as measured at 
launch time. To produce the cleaned and calibrated event files, all the data were reduced 
using the Swift xrtpipeline task (version 11.5) and calibration files from the CALDB 2.8 
releaseQ A cleaned event list was generated using the default pipeline, which removes the 
effects of hot pixels and the bright Earth. The initial data recording was in Windowed Timing 
(WT) mode, due to the large count rate (initially ^1000 count s~^); from approximately 6160 
s after the trigger to the end of the observation (roughly 40ks after the BAT trigger) data 
were taken in Photon Counting ( PC) mode. The X -ray light curve was produced using the 
Swift-XBT light curve generator (JEvans et al.l 120071 ). Note that this tool takes into account 
the pile-up correction which is applied to the data when necessary. From the cleaned event 
list, the source and background spectra were extracted using xselect. We only use grade 0-2 
events in WT mode and grade 0-12 events in PC mode, which optimize the effective area 
and hence the number of collected counts. The ancillary response files (ARF) for the PC and 
WT modes were produced usi ng xrtmkarf, taking into account the point-spread function 
correction ( jMoretti et al.l 120051 ) as well as the exposure map correction. The exposure maps 
were produced using the xrtexpomap task (version 2.4). The vOlO response matrix files 
(RMFs) (swxwTOT02S6_2ooioioivoio.RMF in WT mode and swxpcotoi2S6_2ooioioivoio.rmf in PC 



mode) were used to perform the spectral analysis (JGodet et al.ll2009[ ). 



During the process of creation of the spectra and the light curve, we took great care to 
correct the data for the effect of pile-up in both PC and WT modes. Indeed, the effects of 
pile-up result in an apparent loss of flux, particularly from the center of the point spread 
function (PSF), and a migration of 0-12 grades (considered as grades resulting from X-ray 
events) to higher grades and energies at high count rates. The high observed count rates 
in WT mode (several hundred counts s~^) and PC mode (above 2 counts s~^) indicate that 
the data are piled-up. In WT mode, we produced a sample of grade ratio distribution using 
background-subtracted source event lists created by increasing the size of the inner radius of 
an annular extraction region with an outer radius of 30 pixels (1 pixel = 2.36"). The grade 
ratio distribution for grade 0, 1 & 2 events is defined as the ratio of grade 0, 1 & 2 events 
over the sum of grade 0-2 events per energy bin in the 0.3-10 keV energy range. We then 
compared this sample of grade ratio distribution with that obtained using non piled-up WT 
data in order to estimate the number of pixels to exclude, finding that an exclusion of the 



See [http://heasarc.gsfc.iiasa.gov/docs/swift/analysis/ 1 
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innermost 6 pixels was necessary w hen all the WT data a re used. In PC mode, we used the 
same methodology as described in IVaughan et al.l (120061 ) to estimate the number of pixel; 



s 



which needs to be excluded. We fi t the wings of the so urce radial PSF profile with the XRT 



PSF model (a King function; see iMoretti et al.l l2005l ) and a background component (here 



assumed to be a constant). The exclusion of the innermost 5 pixels enables us to mitigate 
the effects of pile- up. 



2.4. Swift/UVOT Data 



The Ultra- Violet/Optical Telescope (UVOT; iRoming et al.l ( 120051 )) began observing EV 
Lac, at UT 05:16:35 (157s after the Swift-BAT trigger, TO), and performed a sequence 
of automatic observations. The first observation was an 8.08s settling exposure, which is 
usually taken in event/image mode with the v filter while the UVOT stabilizes and prepares 
for settled observations. The second observation, a 200s white finding chart exposure, was 
not completed, because the safety circuit tripped. The trip occurred at UT 05:16:51.3, 
1.4s after the UVOT had finished moving the filter wheel to white. The safety circuit 
tripped due to a rate of >300 000 count s~^ being observed over 31 frames (each of 11ms), 
implying that EV Lac must have been brighter than 7th magnitude in white for a minimum 
of 1.4s. The UVOT resumed observations at ~ 11.89 ks, taking multiple exposures using 
the 7 filters of the UVOT {white, v, b, u, uvwl, uvm2, uvw2). Table [T] lists the observed 
magnitudes measured with UVOT at various times during the fiare. The errors are purely 
statistical and do not refiect any systematic or calibration uncertainties. A description 
of the effective wavelengths and bandwidths of the Swift/UVOT filters can be found at 
http : //swift . gsf c . nasa . gov/docs/swif t/about_swif t/uvot_desc . html. 

Emission in the v band from EV Lac was saturated during the settling exposure, pro- 
viding a limit for the brightness in v of <11.48 magnitude. However, during this exposure, 
EV Lac was sufficiently bright that it was accompanied by an out-of-focus halo of light that 
has undergone multiple reflections within the instrument. Because the count rate of the halo 
and the flux of the source are related linearly, we are able to use the halo to determine the 
brightness of EV Lac during the settling exposure. 

In order to provide a calibration for the halo of EV Lac, we searched the Swift archival 
for stars of known brightness, which are bright enough to produce a halo. We identifled 
9 stars with UVOT observations in the v fllter, with v magnitudes in the GSC-2 database 
ranging from 6.07 to 10.51, and with suitable halos. 



^http://heasarc. gsfc.nasa.gov/cgi-bin/W3Browse/swift. pi 



To determine the count rates of the halos in the UVOT images, source regions and 
background regions were created for EV Lac and for each of the 9 stars. For each object, 
the source region was chosen to cover the entire halo, but excluding the central object itself 
and its diffraction spikes, any visible background sources, and the diffraction troughs in the 
halo. The corresponding background region was taken as a set of circular apertures, placed 
on source-free areas of the sky, with a total area of at least twice that of the source region. 

For each of the 9 stars, the count rate of the halo was determined and then scaled to the 
source area of EV Lac. These values were then plotted against the count rates that would 
be predi cted in the absence of coincidence loss in 5" radius apertures using the v zero point 
given in iPoole et al.l (|2008[ ). A linear function was fitted to obtain the ratio between the 
halo count rates and predicted count rates of the stars. We multiplied the halo count rate of 
EV Lac by this ratio to obtain the corresponding coincidence-loss-free count rate of EV Lac 
in an aperture of 5^^ radius, and converted this into a v magnitude using the zero point of 



Poole et al.l (120081 ). In order to determine the calibration uncertainty on the v magnitude, 
we divided the count rates of the halo calibration stars by the corresponding values of the 
best-fit linear model and then took the standard deviation of the resulting numbers to be 
the la fractional error. This leads to an estimate of f max =7.2±0.2, compared to a V of 10.1 
during quiescence. 

For the rest of the images, in order to be compatible with the UVOT calibration 



( IPoole et al.ll2008[ ). source counts were obtained using an aperture of radius of 5" and back- 
ground counts were obtained from a region of radius 20", positioned on a source-free area of 
the sky. EV Lac was found to be saturated in the white, v and b filters. The software used 
to determine the source count rates can be found in the software release HEADAS 6.5 and 
version 20071106 of the UVOT calibration files. Count rates were corrected for coincidence 



loss and converted into magnitudes following the prescription in iPoole et al.l ( l2008l ) . 

We used the Swift UVOT tool U to estimate the count rate and therefore magnitude of 
the star during quiescence (for an M5 spectral type) in the v, white, and u filters, in order 
to determine the flare enhancements of the star. We obtained additional ToO observations 
of EV Lac nearly a year later to determine its quiescent magnitude in the uvwl, uvw2, and 
uvm2 filters. These quiescent magnitudes are listed in Tabl^ and were subtracted from 
the flare magnitudes to investigate the photometric variations at late times in the flare, as 
illustrated in the bottom panel of Figure [H 



^http://www. mssl.ucl.ac.uk/www_astro/uvot/uvot_observing/uvot_tool. html 
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2.5. Liverpool Telescope 



The Liverpool Telescope took a series of 3x10 s images in the SDSS-r filter ~ 15 minutes 
after the trigger as part of the standard GRB follow-up. No further observations were taken 
once the trigger was identified as a non-GRB. Photometry was done with GAIA calibrating 
SDSS-r images against the USN0B1-R2 magnitudes of the few stars in the field. The night 
was not photometric so systematic errors of 0.2 mag (1 a) are added to all measurements, 
which also account s for the difference between the R2 magnitude of the USNOBl catalog 
(JMonet et al.l 120031 ) and the r' band. Photometry was performed on each of the 3 frames 
and a co-added 30 s frame using an aperture with a radius of 8 arcsec. The differences 
between magnitudes in the three 10 s frames were <0.05 mag so we decided to estimate the 
final magnitudes on the co-added frame. This results in the R2 magnitude of EV Lac being 
7.4±0.2. There is a nearby star, Gl 873B, which at the time of the observations was ~0.5 
arcminutes from EV Lac, which was used to check the photometry. Its magnitude in the 
R2 band is measured to be 10.9±0.2, and its catalogued magnitude is R2~10.6, consistent 
within the errors with the value obtained here. The estimated magnitude for EV Lac during 
this observation is ~ 2 magnitudes brighter than the catalogued one (R2~9.5). 



3. Analysis 



3.1. BAT Spectral Analysis 



We have performed a spectral analysis using XSPEC vl2.4.0 (jArnaudlll996l ) of the BAT 
data for several time intervals during which the hard X-ray spectra were the only data 
collected: viz., during the first slew, first pointing, second slew. During the second pointing, 
(from ^ T0+166:T0-|-961 s) XRT data were also collected. The time intervals during which 
BAT spectra were accumulated are indicated in Figure [1] (b) with blue hori zontal lines. 



Spect ral fitting of the 14-100 keV band used a single temperature APEC model (iSmith et al. 
20011) with abundances fixed to a multiple of the solar values of lAnders fc Grevessd (119891 ) . 



The lAnders fc Grevessd (119891 ) abundance table was used here for bet ter comparison of 



our re sults to other X-ray sp ectroscopic observations of EV Lac, e.g. iLaming fc Hwang 
J2009[ ) and lOsten et al.l J2005h . Because the standard APEC model included in the XSPEC 
distribution only considers emission up to photon energies of 50 keV and the BAT energy 
range extends to 200 keV, we made use of a custom APEC model calculated out to a 
maximum energy bin of 100 keV in order to model correctly the bremsstrahlung contribution 
at high energies. Based on the results from XRT fitting described below, we fixed the 
abundance in the APEC model used to fit the BAT data to Z = 0.4 in order to facilitate 
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comparison between BAT-only and XRT-only spectral fits, the BAT data only being able 
to constrain the product of metallicity and volume emission measure. Figure [2] illustrates 
the variation in plasma temperature, fiux, and volume emission measure (VEM) derived 
from the BAT spectra, and Table [2] lists the plasma parameters and errors derived from 
fitting BAT spectra. The inferred plasma temperature during the first ~100 seconds of BAT 
observations is ~12 keV. To estimate the amount of soft X-ray flux during the times when 
BAT data were collected, we took the best-flt APEC model from flts to the BAT data and 
extrapolated to determine the 0.3-10 keV flux. These are indicated in open red circles in the 
upper right panel of Figure |2j The errors shown are 90% confldence intervals, equivalent to 
±1.6ct. 



3.2. XRT Spectral Analysis 

3.2.1. Time Variation of Plasma Parameters 

In order to investigate spectral evolution during the X-ray flare, we performed a detailed 
time-resolved analysis in WT mode. The statistics in PC mode after pile-up correction are 
not good enough to allow us to perform a meaningful time-resolved analysis. We divided the 
WT data into ten time bins so that we generated in each of the time bins a high statistical 
quality source spectrum with 2 x 10"^ counts after correcting for pile-up. The length of 
the time bins is variable, ranging from 136-521 seconds. Note that the inner radius of the 
annular extraction region used to create the spectra decreases with time since the count rate 
also decreases and the effects of pile-up become less and less dominant. All XRT spectra 
were binned to contain more than 20 counts bin~^, and were flt from 0.3 to 10 keV within 



XSPEC vl2.4.0. ( lArnaudI 119961 ). We performed a single temperature fltting of the data 



using the APEC scaled solar abundance model, as well as spectral flts using two dis crete 



temperature (2T) APEC components. The abundances of lAnders fc Grevessd (119891 ) are 
used here as well. We neglected the absorption column, since the real column towards this 
very nearby star is probably too low to be very well constrained in the 0.3-10 keV band. The 
time variations of the derived parameters (plasma temperature. X-ray flux, volume emission 
measure, and global metal abundance) for the 2T flts are shown in Figure [2J Table E] lists 
the plasma parameters and their errors derived from fltting the XRT spectra. Although the 
errors on the derived abundances during the time slices are large, there is no evidence for a 
signiflcant change in the global metallicity during the flare decay; for the present purposes 
we take the abundance Z to be ~ 0.4. The quiescent X-ray flu x (0.5-7 keV) from Chandra 



observations is ~ 5.8x10 ^^ erg cm ^ s ^ (lOsten et al.ll2005l ). and so the peak measured 



XRT flux is ~ 4200 times this, and even when the flare has declined in flux during the last 
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time bin obtained in PC mode, the flux inferred is still a factor of 10 larger than the Chandra 
value. 



3.2.2. The Emission Line at 6.4 keV 

The 2T APEC spectral fits to the WT mode XRT spectrum obtained from T0+171:T0+279^ 
s revealed positive residuals around 6-7 keV which look like a line feature not included in 
the APEC line list. To investigate that, we added a Gaussian compon ent to the 2-T APEC 



model. The resolution of the XRT at 6 keV was ~ 140 eV at launch ( lOsborne et al.ll2005l ). 
but the CCD has experienced degraded resolution since then due mainly to radiation dam- 
age. The RMFs for WT mode data do not contain proper account of this line broadening 
at high energy due to lack of calibration information. The width of the Gaussian was fixed 
to a value of a = 0.075 keV based on evolution in PC mode dat a of the line width from 
the onboard calibration ^^Fe 5.9 keV line source (see discussion in lGodet et al.l 120091 1. The 
addition of the Gaussian results in a fit improvement of Ax^ = 21 for two parameters and 
does not change the other fitting parameters. The line centroid {E = 6.42 ± 0.04 keV) may 
be consistent with a fluorescence Fe Ka line. Note that the Fe Ka fluorescence line is really 
a doublet with energy separation ~12 eV, but this energy difference is small compared with 
the energy resolution at the time of the observations, and so our use of a single Gaussian 
with a width of 75 eV, broadened by the instrumental resolution of ~ 140 eV is sufficient. 
The use of an F-test gives a probability of 8.5 x 10~^ suggesting that the line is real. How- 



ever, 



Protassov et al.l (120021 ) have opined that the F-test should not be used to assess the 



significance of a line. 

To correctly work out the line significance, we used instead Monte-Carlo simulations 
based on codes provided by C. Hurkett using the posterior pr edictive ?3-value meth od. A 



detailed description of this method can be found in Section 3.2 in lHurkett et al.l (|2008[ ). Here, 
we describe briefly how the codes work. The steps are: i) compute the spectral parameters of 
the model used to fit the data; ii) use these parameters with their errors to generate A^ sets 
of random parameter values and each of these is used to generate a simulated spectrum; iii) 
fit the A^ simulated spectra with the model without a Gaussian line and create a table of the 
X^/d.o.f. values; iv) re-do the same as in iii) but including a Gaussian line, the energy centroid 
varying around 6.4 keV (assuming that the feature is due to a fiuorescence Fe Ka line) and 
the width being fixed at the spectral resolution at the time of the observation {a = 0.075 
keV). A new table of the x^/d.o.f. values is created; v) compute the Ax^ distribution from 
the A^ simulations and compare it with the Ax^ value obtained from the observed spectrum. 
The line significance is then given by the number of times that the simulations produced 
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a Ax^ value larger than the observed one. In order to work, this method needs to run a 
large number of simulations (here, we used A^ = 10^). We used the F-test probability to 
identify which time intervals should be explored further using this method. When the F-test 
probability is low, the above method will not return a high probability, but is needed in 
situations where the F-test gives a high probability; then the posterior predictive p-value 
method can determine whether the line is real. Table H] lists the equivalent widths and fluxes 
of the 6.4 keV feature, along with the results of the F-test and Monte Carlo simulations. We 
performed Monte Carlo simulations for the first six time intervals; the F-test probability for 
time bins 7,8, and 10 were low so we did not do the additional Monte Carlo simulations for 
these intervals. Time bin 9 (T0-(-2215:T0-|-2517 s) shows the presence of a line with an F-test 
probability of 98.5%, so we performed Monte Carlo simulations for this time interval. The 
probability from Monte Carlo simulations is only 97.3%, which we do not consider significant 
(we only considered probabilities > 99% to be statistically significant). Figure [3] shows the 
Ka region of the XRT spectrum for the first five time intervals of the flare decay. 

The fit of the PC spectrum from ~ T0+6.2xl0^:~ T0+ 4.3x10'' s using a two tem- 
perature APEC model gives a relatively good fit (x^/d.o.f. = 236.3/175). The addition of 
a Gaussian line as in the WT spectrum does not improve the fit, not surprising given the 
much lower flux in this dataset. This indicates that the Ka line is not present during this 
late stage of the flare decay. The variability of the Ka line strength is discussed further in 
§3.6. 

The iron Ka line is detected in a variety of X-ray emitting stars, from the solar 
corona to conipact stellar objects to actively accreting pre-main sequence stars (lBailll979 



P 



Giidel fc Nazell2009l ). The formation is generally attributed to a fluorescent process, in which 
an X-ray continuum source shines on a source of neutral or low ionization states of iron (Fe I- 
Fe XII), photoionizing an inner K shell electron and causing the de-excitation of an electron 
from a higher level at this energy. The utility of this line stems from the dependence of the 
line strength on the solid angle subtended by the source of neutral or near-neutral iron. In 
the stellar context where the X-ray continuum emission arises from a loop-top source, and 
the fluorescing material is photospheric iron, the detection of this line thus allows for an 
independent constraint of the height of the flaring loop above the photosphere. The main 
contributors to the observed flux in the Ka line are the total photon luminosity above the 
Ka ionization threshold of 7.11 keV, the fluorescent efficiency of this radiation, and the 
photospheric iron abundance. The fluorescent efficiency changes with photon energy due 
to the energy dependence of the absorption cross section, and so the spectral distribution 
of photons, characterized by the plasma temperature T of the X-ray-emitting material, is 
important for describing the fluorescence process. The fluorescent efficiency also depends on 
the height of the X-ray continuum source above the fluorescent source (this arises from the 
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solid angle dependen ce). The observed flux also d epends on the astrocentric angle (termed 
heliocentric angle in iDrake et al.l l2008t lBailll979l ): the flux is a maximum at zero, which 
corresponds to the center of the stellar disk, and is a minimum when the source is located 
on the stellar limb (astrocentric angle of 90°). The observed flux is written as 



Ka 



f{9)T{T,h)Nr,u 
And? 



photons cm "^ s ^ 



'V. 



where iV7.11 is the luminosity above 7.11 keV, f{9) is a function which describes the angular 
dependence of the emitted flux on astrocentric angle, F is the fluor e scent efficiency, and 
d is the distance to the o bject. This is equ ation 4 in iDrake et al.l (120081 ). We use the 
values in Tables 2 and 3 of iDrake et al.l (l2008[ l for the depende nce of F on t e mpera ture and 
height, and a parameterization of f{0) for different heights Q iDrake et al.l ( 20081 ) perform 
calculations using an iron abundance of Fe/H =3.16x10"^. As we do not have a constraint 
on EV Lac's photospheric iron abundance, we use this default value. We use the best fit 
APEC model parameters for each time interval and create a dummy response from 7-50 
keV to compute the photon density spectrum above 7.11 keV, Nth/ {At: d'^). We use the 
F calculated for a plasma temperature of 100 MK (logT = 8.0), as this is closest to the 
temperatures determined from spectral fltting. Figure H] depicts the observed Ka fluxes and 
the modeling using a point source illumination of the photosphere; the curves falling within 
the observed range of Ka flux provide a constraint on the loop height during the flare, and 
are generally consistent with a loop height less than k. a few tenths of the stellar radius. 
Signiflcant variation of the Ka flux in these time intervals was also found; this is discussed 
in §111 



3.3. XRT+BAT spectral analysis 

During the second pointing, both XRT and BAT data were accumulated, and a joint 
spectral flt was done. For the BAT data this time interval is T0+165.5:T0+961 s, while 
for the XRT this time interval is T0+171:T0+961 s. We added a third temperature APEC 
component to explore whether there was any evidence for other temperature plasma. The 
addition of a third temperature to the model flt reduced residuals below energies of 1 keV, 
and converged on a temperature less than the original two temperatures. Table E] displays 
the flt results from 2- and 3-T APEC models to both XRT and XRT+BAT spectra. Figure E] 
displays the 3-T APEC model flt to the XRT and BAT spectra. No subsequent improve- 
ment in the flt was gained by adding another spectral model, such as a higher temperature 



*There is a typo in their Table 3 for /i = 0.5R* - the as coefficient should be —6.63 x 10^ 
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bremsstrahlung component or a nonthermal hard X-ray emissi on component, as was found 



for the large II Peg flare detected by Swift (jOsten et al.l 120071 ) . A comparison of the fit to 



only the BAT spectrum reveals a higher temperature (10.0^2 s keV, 90% error range) than 
the larger temperature in the two-temperature fit to the XRT+BAT data (6.661q;24 keV, 90% 



error range) although the 90% error ranges of these are only slightly inconsistent. The high 
temperature returned from the XRT-I-BAT spectral fit is the same as that returned from a 
fit to only the XRT data. Thus while the highest temperature measured in the flare during 
the coverage by the Swift satellite appears to be ~ 12 keV (or in units of lO^K, Tg =140), 
the systematic effect of only having BAT data to constrain the temperature may skew this 
value high. 



3.4. UV and Optical Data 

The UVOT and Liverpool Telescope data together constrain the time-scale of the flare 
at optical wavelengths - the Liverpool Telescope data shows that the flare lasted at least 
15 minutes, as the R fllter observations showed a still elevated flux at this time. At late 
times in the flare decay there is only a small (< 0.5 magnitude) offset above the quiescent 
magnitudes, with the u, uvwl, and uvm2 fllters indicating a transient increase at around 
TO+31000 s. The latter is largest in the uvm2 filter, showing an enhancement of 1.2 mag, 
and the increase is restricted to filters sensitive to wavelengths longer than about 2200 A. 

Because the white filter has a response in the ultraviolet, and white-light stellar fiares 
have most of their enhancement in this wavelength range, we use the upper limit obtained in 
this filter to provide a constraint on the fiare area. Previous detailed studies of white light 
stellar fiares have determined that a black-body function with a temperature of ~ lO'^K often 
provides a good fit to the observed broadband fiux distribution, with an area covering factor 
X significantly less than unity. The observed UV/optical fiux contribution due to such fiares 
can be written as 

F, = -^nB,{TBB) (2) 

where Fx is the measured fiux at effective wavelength A, i?^ the stellar radius, d the distance 
to the star, Tbb the black-body temperature, and -Ba(7bb) the Planck function evaluated 
at wavelength A and temperature Tbb- We used the WebPIMMq^l tool to calculate the 
integrated flux in the white fllter, with an input black-body temperature of lO^K, using 



'''http://heasarc.gsfc.nasa.gov/Tools/w3pinims.htnil 
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the count rate of 300,000 counts s ^ at which the UVOT detector safed. The flux in the 
1600-8000 A wavelength range from T0+172.9:T0+174.3 s is >5.5xl0-^ erg cm'^ s'K We 
estimate the average flux over this wavelength range to be >8.6xl0~^^ erg cm~^ s^^ A~^, 
and an area covering factor x of > 0.03 for the effective wavelength of the white filter 
3471 A. For the radius of EV Lac, R^, is roughly 0.36 R© (see discussion in §3.5), and this 
works out to an area of >2xl0^^ cm^. We also calculated the range in x for black-body 
temperatures from 800 K to 20,000 K, find ing that x varies from > 0.08 to 0.003 for this 



range of temperature. iHawley et al.l ( 120031 ) determined Tbb ~ 9000 K and x ~ 10~^ for 
several flares on the nearby flare star AD Leo, finding that each flare in their study had a 
similar blackbody temperature regardless of flare energy, peak flux, or total duration. If we 
assume that this flare had a similar blackbody temperature, it indicates an area covering 
factor ?^480 times the well-studied, but smaller amplitude, flares seen on AD Leo. 

We used the WebPIMMS tool to calculate also the integrated flux in the v filter at 
T0+156.58:T0-|-164.78 s with an input black-body temperature of lO^K, and converting 
the observed v magnitude of 7.2±0.2 into a count rate using the zero-point of 17.89. The 
5000-6000 A flux from T0+156.58:T0+164.78 s is 5.70 xlO'^ erg cm'^ s-\ while the 0.3-10 
keV X-ray flux near this time (T0+171:T0+306 s) is 2.45x10"^ erg cm~^ s~^, leading to 
an X-ray /optical flux ratio fx/fv of 4x10^'^. The quiescent X-ray flux is about 5.8x10"^^ 
erg cm~^ s~^, and the flare enhancement was 2.9 magnitudes in the v filter, allowing us to 
compare this ratio with quiescent values. It is about 280 times higher than in quiescence. 
These values may assist others in trying to make identifications of likely flare star targets 
based on optical and X-ray flux ratios. 



3.5. Hydrodynamical modeling of the flare decay 



Reale et al.l (119971 ) have described a method to estimate the coronal loop length from 
the decay of a stellar flare. The technique uses temporal information on the rate of decrease 
of the flare intensity during the decay phase, as well as time-resolved spectral analysis to 
describe the time dependence of the flaring plasma's temperature and density (which can be 
estimated from the time dependence of the volume emission measure, under the assumption 
that the flare geometry and hence volume does not change appreciably). Most flares exhibit 
an exponential decay, where the decay time is determined by the balance between cooling 
from radiation and conduction, as well as the presence or absence of sustained heating. 
The slope of data points in the temperature-density plane during the flare decay traces the 
heating time scale, which can be used to estimate the loop length for the assumption of a 
single semi-circular loop with constant cross section. By performing hydrodynamic modeling 
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of decaying flaring loops with a range of dimensions and time-scales, and folding the modeled 
flare emission through the response of a specific instrument, the method can be applied to 
the particular bandp ass and instrumen tal response. This approach has been calibrated with 



resolved solar flares (JReale et al.lll997l ) showing good agreement. It has also been applied to 
a few singul ar cases of stellar flares with additional leng th scale constraints, showing general 



consistency (JFavata fc Schmittlll999l : iTesta et al 



n ai ieng t 

D boosh . 



This method has been applied to the Swift /XRT bandpass from 0.3-9.5 keV for a global 
coronal abun dance Z of 0.4 tim es the solar value. The thermodynamic loop decay time can 



be expressed ( ISerio et al.lll99l[ ) as 



Tth = "VV ^max (3) 

where a=3.7xl0~^ cm~^ s^^K^/^, / is the loop half-length in cm, and T^ax is the flare max- 
imum temperature (K), 

T^ax = 0.0261T,^br (4) 

and Tobs is the maximum best-flt temperature derived from single temperature fltting of the 
data. The ratio of the observed exponential light curve decay time tlc to the thermody- 
namic decay time Tfh can be written as a function which depends on the slope ( of the decay 
in the log(ne-Te) plane (or equivalently, \og{VVEM — T) plane): 

rLc/rtH = -^ + qa = F{C) • (5) 

The parameters ga, Ca, and Ca are determined for the Swift/XRT instrument to be ga=0.67±0.33, 
Ca=1.81±0.21, and Ca=0.1±0.05. Combining the above expression with the one for the ther- 
modynamic loop decay time, a relationship between flare maximum temperature, light curve 
exponential decay time, and slope in the density-temperature plane can be used to estimate 
the flare half-length: 



7 'LC V -^ max / « \ 

valid for 0.4 < C ^ 1-9. The uncertainty of the loop length is determined by propagating 
errors from the tlci C; ^obs values determined from data, and the uncertainties of Ca, Ca, and 
qa from the numerical calibration of the method to Swift/XRT data. 

Section 13.21 discussed the time-resolved spectral analysis of the XRT data. In order to 
obtain a density-temperature path, the method requires single values of flare temperature 
and emission measure which can be obtained with a IT fltting. The same fltting is applied 
to the flaring loop simulations used to calibrate the method. In this case, the flt provides 
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an effective temperature weighted average on the emission measure components. The slope 
( in the log{VVEM — T) plane using the isothermal fits is C = 0.74±0.03 (as shown in 
Figure E]), with a maximum observed temperature Tobs of 67±4 MK. This corresponds to 
^max of 142 MK. Flares in which the time-scale for the decay of plasma heating is gradual 
have small values of (, near 0.4, whereas flares in which the heating is shut off abruptly have 
larger values of (, near 1.9. The exponential decay time of the XRT light curve for the first 
~2800 s is tlc =1014±5 s. These can be combined using the equations above to arrive at a 
loop semi-length estimate / of 9.3±1.4xl0^ cm. The relatively short flare duration and the 
slope C signiflcantly larger than the low asymptot ic value (0.4 ) both support a single loop 
flare rather than a multiple-loop two-ribbon flare ( lRealell2003l ). Figure |6] displays a steeper 
value of ( in the initial stages of the flare decay, closer to a value of 1.1. Discussion in §4.1 
showed that the flare decay had an initially fast decline in the flrst 1000 s, followed by a 
slower decline (r^c of 647 s followed by tlc of 1042 s). Using these values of ( and r to 
determine the loop semi-length gives consistent results with the method above. 

Expressing the loop length estimate as a fraction of EV Lac's stellar radius requires 
a precise estimate of the star's radius. Based on m ass-luminosity rela. tionship arguments, 
EV Lac's mass is likely ~ 0.35 Mq (see discussion in iFavata et al.ll2000l ). EV Lac's spectral 
type of M3.5 places it at the dividing line separating fully convective stars from those with 
a thin inner radiative zone. The exact radius will depend on metallicity and possibly also 
on other factors: recent work has shown a systematic discrepancy of stellar radii of low 
mass stars determined from eclipsing binaries compared with l ow-mass stellar isochrones 
( iFernandez et al.l l2009l . and references therein). The models of IChabrier fc Baraffd (Il997| ) 
give a stellar radius of 0.36 R© for a 0.35 M© solar metallicity star older than about 300 MYr. 
We add in an uncertainty of 10% to this value of stellar radius to account for any systematic 
effect not included in the model determination. With this, the loop semi-length estimate 
expressed as a relative length is //-R* = 0.37±0.07. If the flaring loop has a semicircular 
geometry and stands vertically on the stellar surface, then the maximum loop height h can 
be expressed as /i/i?*=0.24±0.04. Since the loop might be inclined on the surface, this value 
of h should be considered as an upper limit. 



3.6. Fe Ka 6.4 keV and Its Variations 



To date there have been only a handful of detect ions of Fe Ka 6.4 keV emission in 
stars without disks ( lOsten et al.ll2007l : iTesta et al.ll2008l ). and these have occurred in evolved 
stars. The best studied case, due to our abihty to spatially resolve flaring regions, is the Sun 
( lBailll979l ). where the hot continuum X-ray radiation from a flaring coronal loop-top source 
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illuminates the photosphere and causes fluorescence of photospheric iron. In addition, there 
have been more detections of Fe Ka emission in young stars, where the possible presence 
of a disk comphcates the geometry, with reflection o f the stellar X-ray con tinuum off the 



surrounding disk likely bei ng a dominant contri butor (JTsujimoto et al.ll2005l ). EV Lac does 



not appear to have a disk (JLestrade et al.ll2006l ). and thus the detection of the fluorescence 



line is the first time this has been seen in a main sequence star other than the Sun. 

The Ka feature is detected most significantly in the early stages of the flare decay, and 
so we investigated the time-dependent characteristics and formation of the line. Due to the 
need to mitigate the pile-up effects discussed in §2.3, there is not enough signal on shorter 
time intervals than the time bins we have used in Table HI Extracting spectra over larger 
time intervals does not improve the significance of Ka equivalent width much, indicating 
that it is likely fluctuating on shorter time-scales compared to the size of the larger bins. In 
the first five time intervals the photoionizing hot continuum, defined as A^7.ii/47r(i^, declines 
by a factor of 6, while the plasma cools from 8.4 keV to 4.8 keV. This results in a general 
decline in the strength of the Ka feature. The Ka line fluxes through the decay phase of the 
flare are generally consistent with loop heights on the order of a few tenths of a stellar radius. 
The hydrodynamic modelling of § 13.51 returned an estimate of the height of the flaring loop 
of h/R^ =0.24±0.04. Thus the two methods are in agreement on the coronal loop height 
within the uncertainties of factors of a few. 

The time variation of the Ka line flux is generally consistent with fluorescent emission 
from compact loops and the cooling of the flaring plasma occurring during the decay of the 
flare. However, the initial stages of the flare decay show some marked variations on short 
time scales. For 135 s (T0+171:T0+306 s) the Ka flux is detected at >99% confidence, 
indicating a compact loop of height <0.1 R^, and then in the next 202 s (T0+306:T0-|-508 
s) the line is undetected, with a lower limit of zero flux. In the third time segment, spanning 
163 s (T0+508:T0-(-671 s), the line flux is again present but at a level greater than what 
is predicted based on fluorescence from the derived loop height — in fact, the observed 
Ka flux is not consistent with any set of fluorescent models, even for those that go to zero 
height. In this third time interval, the most compact loops (with h asymptoting to zero 
height) would have a Ka line fluorescent flux which is ^8xl0~^ photons cm~^ s~^. The 
strength of the Ka flux is about three times larger than that, with an total flux of about 
2.64x10^^ photons cm~^ s~^. The excess emission is conservatively in the range 0.0089 < 
FeAkeV < 0.028 photons cm^^ s^^ (using the 3a error bars in Table Hj). The time-scale for 
variation of this flux is rather short, being present for only 163 s. The duration is short 
enough that changes in astrocentric angle (such as may happen during a large change in 
phase of stellar rotation) can be excluded. The fluorescent effi ciency F depends on the 



fluorescing spectrum, flare height, and relative iron abundance. iDrake et al 



([20081) show 
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that r varies very slowly with temperature, decreasing by a factor of 2 over a factor of 30 in 
temperature. The temperatures at which the Ka feature is detected are 6-8 keV (70-90 MK) 
and the integrated photon density spectrum decreases by a factor of 3 during the first three 
time intervals. The r elative fluorescent efficiency varies with abundance by weaker than a 
proportional relation (JDrake et al.l 120081 ) . and there is no evidence for abundance variations 
taking place during the decay of this flare. 

One possibility that has been used to explain Fe Ka emission during the impulsive phase 
of a few solar flares is the association o f the Kg line flux v ariations with hard X -ray bursts 
jTanakaet al.lll984J : ISmslie et aDll986l : IZarro et al.lll992[ ). iTanaka et all Jl984l ) attributed 
the intense Ka emission in the early phases of the solar flare to fluorescence from the power- 
law hard X-ray bremsstrahlung emission, noting that further Ka emission af ter the decay 



of the hard X-ray flux was consistent with fluorescence from a thermal source. lEmslie et al. 



(119861 ) determined in the solar flare they considered that the excess Ka emission above 
that produced by therma l photoioniz a tion w as consistent with coUisional ionization from 
accelerated electrons, and IZarro et al.l ( 1l992l ) noted during the impulsive phase of the solar 
flare they analyzed that a thermal photoionization interpretation of the observed Ka line flux 
could not account for all of the observed Ka line flux, requiring a component from coUisional 
ionization. In contrast to the situation presented in these solar flares, in the present case we 
are detecting Ka line emission during the decay phase of the stellar flare, and the variability 
appears during this decay phase. As noted in § 13.31 the analysis of the combined XRT and 
BAT spectra during the flare decay over the time interval T0+171:T0+961 s showed no 
evidence for an additional spectral component beyond the two temperature APEC model 
and the Ka line emission at 6.4 keV. The BAT data did not have enough counts to extract 
spectra over shorter time intervals, and so we use this spectra to constrain the action of 
nonthermal electrons over the larger time interval. 



Following the discussion in lEmslie et al.l (119861 ). we can determine the amount of Fe 
Ka emission produced by coUisional ionization. The observed Ka flux depends on the 
characteristics of the distribution of accelerated electrons, namely the power-law index 6, 
the total power in the beam Fq, and the low-energy cutoff Eq. Additionally, the Ka flux 
depends on the column density in the Ka -emitting layer (iV*), as some fraction of accelerated 
electrons will lose their energy via collisions in dense regions. We computed the expected Ka 
flux on a grid on 6, Fq, and A^*, and cutoff values Eq =20 keV and Eq =100 keV, and found 
ranges of parameter space where the Ka flux due to coUisional ionization was at the level of 
the Ka flux observed in the third time interval — the measured Ka flux and errors minus 
the maximum amount expected from a fluorescence process (0.0089 < FQ^keV (photons cm~^ 
s~^) < 0.028). We then used these parameters to estimate the amount of nonthermal thick- 
target bremsstrahlung flux. Since the observed flux level at 20 keV during the time interval 
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T0+171:T0+961 s is ^10^^ photons cm"^ s~^ keV~^ and there is no evidence for nonthermal 
contribution to this flux, we used this value to determine plausible ranges of parameters which 
could both produce enough Ka emission to provide the required extra Ka flux and not be 
detected in the BAT spectrum. For column density in the Ka -emitting layer of 10^^ cm~^ 
or less, low energy cutoff Eq of 100 keV, values of total power in the beam between 10^° and 
]^g33 gj,g g-i^ g^j^^ ^ between 6 and 8, the photon flux at 20 keV is < 10~^ photons cm~^ 
s~^ keV~^. Figure [7] displays the range of parameter space, with pluses indicating where 
the expected flux was within the measurement range for the excess Ka emission line flux, 
and the pluses with circles around them indicate parameter values consistent with a 20 keV 
thick target bremsstrahlung flux of < 10^'^ photons cm~^ s^^. Thus it is plausible that the 
excess emission in the Ka line could be produced by coUisional ionization, and yet not be 
detected as an additional spectral component in the BAT spectr um. There are examp les of 



solar flares whic h have high val ues of the low energy cutoff Eq (jWarmuth et al.ll2009l ) and 



high values of 6 ( ISui et al.ll2004J ). 



The optical area constraints from the white light flare discussed in §3.4 lead to an area 
A > 2x 10^^ cm^ . This is the area of the photosphere involved in the flare, and if we 
interpret this as the foot-point area, then we can combine this measurement with the above 
determined ranges of total electron beam power, to arrive at an estimate of the nonthermal 
beam flux. For a beam power between 10^° and 10^^ erg s~^, this leads to a nonthermal 
beam flux of 10^^-10^'^ erg cm~^ s~^. The values typic ally found in solar flares are 10^° erg 



cm s , with 10 being found in some large flares. lAUred et al.l (|2005f) di scuss radiative 



hydrodynamic models of solar flares, and a later paper (lAUred et al.l 120061 ) extended the 
results to M dwarf stellar flares; they found that for the same beam fluxes, the time-scales 
in both the solar atmosphere and M dwarf atmosphere were shorter for larger beam fluxes. 
For a beam flux of lO^^erg cm~^ s~^, the dynamics of the impulsive phase of the flare lasts < 
15 seconds in both cases. These short time-scales are consistent with the short time-scales 
observed in the flare on EV Lac. 



4. Discussion 

4.1. Flare Energy Estimates 

Because of the heterogeneous wavelength and temporal coverage during the flare decay, 
we made flare energy estimates a few different ways. The Swift observations cover the initial 
~3000 s of the flare decay and additional observations extending up to 4x10^ s after the 
trigger. Table E] lists the determination of radiated energy at various times and wavelength 
ranges during the flare. We estimated the energy radiated in the 0.3-10 keV range during 
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the time intervals in which only BAT spectra were obtained by extrapolating the best fit 
to the 15-100 keV X-ray spectrum into this energy range. We also examined the temporal 
evolution of the 0.3-10 keV flux as shown in Figure [2J Fitting this with a double exponential 
function with a break around 1000 s after the trigger, we find an initial exponential decay 
Ti of 647 s and second decay T2 of ~ 1042 s. We then extrapolate to the point where the 
flux would equal the measured quiescent X-ray flux, ?a5.8xl0~^^ erg cm~^ s~^, which is ~ 
8300 s after TO. In the 0.3-10 keV band we estimate the radiated energy to be ?^7.3xl0'^^ 
ergs from T0:T0-|-8300 s. The flux measured during the time span 6160-19792 s after TO 
is 10 times higher than the quiescent flux, even though extrapolation of the flux due to an 
exponential flare decay would indicate a return to quiescent conditions at around TO-l-8300 
s. This suggests that the flare likely was more complex, with a re-brightening during this gap 
in time, and a continuation beyond the Swift monitoring time. Since this was apparently a 
long-duration flare, our estimates of the radiated flare energy are actually lower limits to the 
total energy radiated in X-rays. The flare radiated energies are at the high end of what has 
been observed for flares on other M dwarfs. Note that the exponential time-scales derived 
from BAT and Konus in equivalent time ranges are shorter: from T0:~T0-|-300 s r is 175^^4 
s for BAT (14-30.1 keV) and 187±7 s for Konus (18-70 keV). 

We note that in the first spectrum of the flare, from (T0-39:T0 s), ~62% of the stellar 
bolometric flux was recorded at hard X-ray energies (14-100 keV). By extrapolating the 
temperature and emission measure determined from that spectral analysis to smaller X-ray 
energies, we estimate that in the XRT band (0.3-10 keV) the flux may have been as high as 
4.1x10"^ erg cm~^ s^^, or ^ 2.4 times the bolomet ric flux from t he sta r. This is about 7000 



times larger than the quiescent flux estimated by lOsten et al.l (120051 ). Thus from 0.3-100 
keV the peak luminosity was 1.6x10^^ erg s~^, or a value Lx/Lboi of ~3.1, from T0-39:T0 
s. The upper right panel of Figure [2] displays the bolometric flux, along with the temporal 
evolution of the flux in the 0.3-10 and 15-100 keV energy ranges; for the first ~ 500 s after 
the trigger the X-ray flux exceeded the bolometric flux. 

At first glance it may seem problematic to have a flare whose luminosity exceeds the 
normal stellar bolometric luminosity. The free energy which powers the flare is built up over 
long time-scales in a more or less continuous manner. However, that energy can be released 
suddenly through explosive magnetic reconnection. The largest e nergy flares which can be 



produced relate to the largest active region size which can exist (jAschwanded 120071 ) . This 
flare was an extreme example of instantaneous magnetic energy release, but the radiated 
energy in the X-ray band is typical of longer duration flares. 
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4.2. Comparison to Other Superflares 



While there have been stellar flares detected above 10 keV in the past (mostly with 
the BeppoSAX satellite), these S wift /BAT flares ( the one on EV Lac reported here and the 
superflare on II Peg reported by lOsten et al.l 120071 ) are different in that they were observed 
using autonomous triggers based on their hard X-ray flux rather than in pointed observations 
which serendipitously detected large stellar flares |j. The two triggered flares differ in several 
important respects: the maximum plasma temperature achieved was very large in the early 
stages of the flare on II Peg (peaking at ~180 MK), and the temperatures measured with both 
XRT and BAT spectra at about 3 minutes past the triggers indicate hotter temperatures for 
the case of II Peg versus EV Lac (120 MK and 100 MK, respectively). The II Peg flare had a 
longer rise time (1000-2000 s) compared with <120 s for EV Lac, and while neither flare was 
observed in its entirety, the decay of the flare on II Peg had a slower decline than the flare 
on EV Lac. Some of these differences are likely due to the different stellar environments: 
II Peg is a tidally locked binary system consisting of a K subgiant and M dwarf companion 
with a likely age of a few Gyr, while EV Lac is a younger single M dwarf flare star, with a 
much higher gravity. 

At its peak, the EV Lac flare is ~2.5 times II Peg's maximum 14-40 keV hard X-ray 
count rate, and the XRT flux is ~ factor of 4 larger. Due to the different distances (42 pc for 
II Peg vs. 5.06 pc for EV Lac), the X-ray luminosity for II Peg during the peak of that flare 
was larger, at 1.3x10^^ erg s~^ (0.8-10 keV), compared with 6.5x10^^ erg s~^ for EV Lac in 
the same energy band. Expressed as a fraction of the star's bolometric luminosity, though, 
the EV Lac flare is larger, with i^o.s-iokev/Lboi ~ 1-2, compared to Lo.s-iokev/Lboi = 0.24 
for II Peg. As discussed in §4.1, up to several times EV Lac's normal bolometric luminosity 
may have been radiated at X-ray energies when estimating the soft X-ray flux in the early 
stages of the flare. Another crucia l difference i s the e xistence of additional emission in the 
hard X-ray band in II Peg, which lOsten et al.l (120071 ) interpreted as nonthermal emission; 
in the EV Lac flare there is no such evidence. Both flares revealed evidence for Fe Ka 6.4 
keV emissio n, but the EV Lac flare sports a larger maximum equivalent width and photon 
flux. While lOsten et al.l ( 120071 ) described the formation of the K g line flux seen in II Peg as 
originating entirely from coUisional ionization, iDrake et al.l ( 120081 ) showed that a fluorescence 
interpretation is more likely. The Ka line flux shows temporal variations in both flares, and 
here we interpret the excess Ka emission as suggestive of the action of nonthermal electrons. 



^ Since this flare, there have been two additional autonomously triggered stellar flares observed with 
Swift: ,http://gcn. gsfc.nasa.gov/gcn3/8371.gcn3, and http://gcn.gsfc.nasa.gov/gcn3/8378.gcn3; these are the 
subject of future papers. 
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Besides the two flares discussed above, two other flares have detecti ons in the 10-50 
keV band made with Beppo SAX — a flare on the RS CVn binary UX Ari (JFranciosini et al. 



2OOII ). and a flare on Algol ( iFavata fc Schniittlll999[ ). The common characteristics of these 
flares were a super- hot temperature (peak temperatures in excess of 100 MK), large X- 
ray luminosity (generally ~10^^ erg s~^, in different energy bands), and long duration, as 
measured by the time to decay from peak count rate to half that value (~ 8 and 12 hours 
for the Algol and UX Ari flare, respectively). The flare on UX Ari showed no abundance 
variations, while the one on Algol did show factor of ~3 abundance enhancements in the 
early stages of the flare. Observations of stellar fl ares with the Gin ga satellite have detected 
X-ray emission out to 20 keV in a few large flares. iPan et al.l (119971 ) found peak temperatures 
near 100 MK and and emission measures near 10^^ cm~^ from a large flare identified with the 
M dwarf star EQ1839. 6-1-8002, and flare peak temperatures fr om 60-80 MK were determined 



for flares on the active binaries Algol, UX Ari, and II Peg (iStern et al.l 119921 : iTsuru et al. 
19891 : IPovle et al.lll99lh . 



ATEL # 1689 and # 1718 @ describe an X-ray transient detected by XMM-Newton 
which may be an M star. With a distance of 400 pc, its peak X-ray luminosity in the 0.2-10 
keV range was 5.5x10^^ erg s~^. The authors determine the spectral type of the star to be 
M2/3, and thus this X-ray flare, if associated with this star, would likely have exceeded its 
bolometric luminosity (~10^^ erg s~^). 



4.3. Comparison to Other Flares on EV Lac 

Prior to the large flare discussed in this paper, Swift/BAT detected lower amplitude hard 
X-ray transient emission from EV Lac on three occasions which were not strong enough to 
trigger an autonomous slew. ATEL #1499 u describes these. The peak signal-to-noise ratio 
in the BAT data of the autonomously triggered stellar flare was 19.47, while two possible 
events on 6 June 2006 and 20 May 2007 had SNR of 3.5 and 3.1 sigma, respectively. An event 
on 22 April 2008 was a 4.0 sigma detection over 64 seconds, and is a likely flare detection. 

EV L ac earlier underwen t a large flare which was detected by the ASCA satellite, dis- 
cussed by iFavata et al.l (I2OOOI ). The maximum plasma temperature diagnosed was 6.3 keV 
{^ 70 MK) with a peak X-ray luminosity (0.1-10 keV) of -lO^i erg s^^ (25% of Lboi), 



^http://www.astronomerstelegram.org/?read=1689 
^http://www.astronomerstelegram.org/?read=1718 
^http://www.astronomerstelegram.org/?read=:1499 
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and radiated energy at X-ray wavelengths (0.1-10 keV) of 1.5x10 erg (^ 300 times the 
quiescent level). Like this flare, the one observed by ASCA was relatively short, with an 
exponential decay of 1800s, and came from a relatively compact flaring loop (/ <0.5R^). 
Visual inspection of the Konus-Wind data for the time of the flare caught by the ASCA 
satelhte shows no signiflcant emission, suggesting the flare detected by Swift was harder 
and more luminous in the hard X-ray range. Another key difference between the flares was 
the presence of signiflcant abundance variations in the ASCA flare (factor of ~6 between 
maximum and minimum Z). 

A mode rate flare on EV Lac see n with the Suzaku X-ray telescope lasting 1500s was 
discussed by iLaming &: Hwangl (|2009[ ). In comparison with the spectrum extracted outside 
of flaring intervals, the abundance during the flare appears to have increased by f actors of 



several. This behavior is similar to what was observed in the large ASCA flare. ISchmitt 



(119941 ) discussed a flare seen on EV Lac in the ROSAT All-Sky Survey which lasted 24 hours, 
with a peak temperature near 25 MK. Numerous small flares have also been observed with 



Chandra — lOsten et al.l (120051) d e scribe d nine small flares, with highest temperatures near 



30 MK, and iHuenemoerder et al.l ( l2009l ) discussed another Chandra observation exhibiting 
numerous low-level flares. 



4.4. Coronal Length Scales 

The determination of coronal length scales provides an important constraint on the 
energetics and dynamics of coronal processes. The use of hydrodynamic modelling to describe 
the variation of plasma temperature and volume emission measure through the flare decay 
phase allows an estimate of the coronal loop length to be made, assuming that a single flaring 
loop dominates the emission. Likewise, the detection of the Fe Ka line allows a constraint to 
be made on the height above the surface of a coronal loop illuminating the cold iron which 
fluoresces to form the Ka line. While the hydrodynamic modelling explicitly assumes that 
a single loop is involved in the flare, the Ka emission line strength depends mainly on the 
height. The fact that the two are consistent to within a factor of a few suggests that the 
former assumption is roughly consistent with the data. 

Loop lengths derived from an analysis of the time variation of flare light curves in con- 
junction with variation of temperature and emission measure using hydrodynamical mod- 
elling of dec aying flaring l oops i ndicates relatively short loops for two large flares observed 
on EV Lac: iFavata et al.l ( 120001 ) derive a loop semi-length I of ~0.5R* during a large flare 
seen with the ASCA satellite, while in the present case a comparable loop length is derived, 
I ~ 0.37R^. 
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Another method to estimate length scales involves resonance scattering in emission lines, 
where the expected flux of density-insensitive resonance lines in emission departs from that 
expected under optically thin conditions. Then the s c atteri ng length can be deduced if the 



electron density in the medium is known. iTesta et al.l (120071 ) determined from time-averaged 



coronal spectra obtained with Chandra's HETGS of EV Lac that the ratio of Lyman-a to 
Lyman-/3 lines of O VII was more than 3a below the theoretical value assuming optically thin 
emission. Using an escape probability formalism they estimated a path length of 1.6x10^ cm 
(0.06R^) for EV Lac's O Vll-emitting material. Such techniques rely upon the measurement 
of individual lines, which is not possible with the XRT on Swift. These length scales are 
smaller than those determined for the flaring loop in this case, but are within an order of 
magnitude. 

If we use the flare area determined from the white light flare and assume that this is 
the area of the foot-point of a single coronal loop, then we can estimate the radius of the 
coronal loop at the base under the assumption of simple cylindrical loop. From §3.4, A 
>2xl0^^ cm~^, so dividing by two for a contribution from each foot-point, we arrive at a 
foot-point radius of r ;=slO^ cm. The loop length estimates from §3.5 give a semi-length of 
/ =9.3x10^ cm, and combining the two gives an aspect ra tio a = r/{2l) of 0.1. This is the 



value commonly assumed in solar flares (JGolub et al.lll980l ). Figure [8] illustrates the relation- 



ship between the X-ray photons, coronal loop, white-light flare emission, and region of Ka 
emission. The plasma volume emission measure (VEM) contains spatial information through 
the total volume, which can be used to constrain loop length given an electron density, or 
used to estimate electron density if a loop length is assumed. With a set equal to 0.1, the 
volume can be written as 

V = SnaH^ (7) 

and the VEM as 

VEM = nlV = 8nl7TaH^ . (8) 

The peak 0.3-10 keV VEM during this flare was ~ 10^^ cm~^. If we use the loop length 
derived from the hydrodynamic modelling to estimate an electron density, and assume that 
there was only one flaring loop, then an upper limit on Ug is around 6x10^^ cm, which is 
quite high. It is an upper limit on Ug because the loop length might be more than 9.3x10^ 
cm and inclined on the surface, so the loop volume might be larger, and hence the density 
smaller. We have no constraint from X-rays on the loop number or cross-section area. If the 
number of loops was larger than one, this would not change the observed plasma evolution. 
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4.5. Astrobiological Implications 



Due to their lower luminosity, M dwarfs have a habitable zone which is 5-10 times closer 
than for a solar-like star. This decreased distance renders an exoplanet more susceptible to 
influences by the parent star. It has already been recognized that flares from dMe stars 
affect the atmospheres of a ny terrestrial exoplanets through frequent stochastic ionizing 



radiation (jSmith et al.ll2004j ). We can quantify the effect that the flare being discussed in 
this paper would have on an exoplanet in the EV Lac's habitable zone by comparing it to 
solar flares using the solar X-ray flare classification. Using the spectral model during the 
trigger observation (T0-39:T0 s), we calculate the flux in the 1-8 A (1.54-12.4 keV) range, 
and scale from a distance of 5.06 pc to 0.1 AU. The solar flare classification is Xn = nx 10~^ 
W m~^, and for this flare the scaled flux is 3600 W m~^, or an X3.6xl0^. For comparison, 
the largest solar flare is about X28. 



Smith et al.l ( 120041 ) also determined that roughly 4% of the ionizing X-ray radiation can 
penetrate a terrestrial atmosphere, through redistribution into biologically damaging UV 
radiation. In the flare considered here on EV Lac, the X-ray emission was approximately 
7000 times the quiescent level during the peak of the flare, suggesting that the UV radiation 
which a putative terrestrial planet around EV Lac would have experienced was increased by 
about a factor of 280 due only to the X-rays. Since the X-ray flare was accompanied by white 
light flare, which was a factor of 76 brighter in the UV than the quiescent UV emission, the 
total increase in the incident radiation on a terrestrial planet would have been ~350 times 
larger. 

The UV flux from a parent star can affect the evolution of a terrestrial planetary atmso- 
phere, and can be bo th damaging and necessary for biogenic processes on a terrestrial planet. 



Buccino et al.l (120071 ) defined a UV habitable zone which balances these two competing pres- 
sures, and found that for inactive M stars the UV habitable zone was closer to the parent 
star than the traditional liquid water habitable zone. During moderate flares, however, the 
temporary increase in UV radiation made the two habitable zones overlap, suggesting that 
there may be a beneficial role for frequent, moderate sized flares in supporting habitability 
on exoplanets around M stars. 

The computation of habitable zone distance is partly determined by the star's bolometric 
flux. The extreme flare here represents a case of a transient increase in the star's bolometric 
flux, with an energy distribution significantly blue-ward of the quiescent spectral energy 
distribution. Although this magnitude of flare is likely rare, it likely represents the "tip of 
the iceberg" of a whole ensemble of more frequent but less energetic flares whose integrated 
emission may change the time-averaged bolometric luminosity over the value estimated from 
the luminosity in the visual band using the standard bolometric correction. This is an 
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additional factor which may influence habitable zone distances for M stars. Future attempts 
to characterize the habitability of exoplanets around M stars should take into account the 
influence of such extreme events. 



Conclusion 



The flare studied in this paper is unique in two main aspects: the level of X-ray radiation 
marks it as one of the most extreme flares yet observed in terms of the enhancement of ~7000 
compared to the usual emission levels, and it is one of only a handful of stellar flares from 
normal stars without disks to exhibit the Fe Ka line, with a maximum equivalent width of 
>200 eV. 

It is somewhat ironic that altho ugh early expla nations for the source of GRBs considered 
gamma ray emission from flare stars (JLi et al.lll995l ). they were quickly swept aside due to the 
extreme nature of the events needed and the apparent discontinuity with the population of 
observed stellar flares. Now that gamma-ray burst detectors have become more sensitive, we 
are indeed finding a small population of very energetic stellar flares which produce transient 
hard X-ray flux at levels consistent with those seen from extragalactic GRBs. These highly 
luminous events do appear to share characteristics with lower amplitude flares, and are likely 
an extension of the same underlying energy release processes. Thus, the fact that these flares 
exhibit Fe Ka line emission and evidence (both direct and indirect) for the existence of 
nonthermal hard X-ray emission does not imply that there is something s pecial about these 
flares; stated in the solar physics literature as the "big flare syndrome" ( lKahlerlll982l ). the 
larger the energy release, the more likely one is to see different flare energy manifestations. 
Perhaps, the most interesting finding is the detection of Fe Ka line emission, as its presence 
in the X-ray spectra of normal stars without disks gives an independent constraint on the 
height of flaring coronal loops, and provides a consistency check with results from flare 
hydrodynamic models. As previous attempts to characteriz e the length of flaring loops 
using different theoretical methods has led to different values (JFavata &: Schmittlll999l ). the 
use of independent methods is an extremely useful 'sanity check' on these models. 

The detection of Fe Ka line emission in flares from normal stars is even more interesting 
in light of its utility in understanding the formation of this line in young stellar objects, where 
the possible presence of a disk complicates the interpretation. In young stellar objects, the 
Fe Ka line is somewhat more commonly detected, but due to the often unknown nature of 
disk geometry arou nd these stars, the int erpretation of the strength of the emission is more 
ambiguous. Indeed, iGiardino et al.l (120071 ) determined that the observed Fe Ka line strength 
in one young stellar object exceeded that allowed either by fluorescence or disk emission, and 
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suggested the presence of nonthermal electrons in the star-disk interaction. The advantage of 
the Swift-triggered stellar flares is not only their high intensity which allows for the detection 
of this line, but the ability to study the line flux variability at high time resolution, along 
with simultaneous constraints on the hard X-ray flux level. 

The size of the flare, in terms of its peak X-ray luminosity exceeding the non-flaring 
stellar bolometric luminosity, is impressive. The existence of such flares can provide impor- 
tant constraints on the time-scales for energy storage and release in a stellar context. The 
multi-wavelength observations which comprise this flare study have enabled the determina- 
tion of both the rise time of the flare (Konus/Wind and Swift/BAT constraints) as well as 
the photospheric contribution to the flare at early times (Swift/UVOT v and white filter 
constraints). While more extreme stellar white- light flares have been observed from other 
M dwarf stars, the combination of the UV, optical, soft and hard X-ray measurements were 
key for determining the possible contribution of nonthermal electron beams to the Fe Ka 
line flux, and confirming the self-consistency of this explanation by estimation of nonthermal 
beam fluxes. While we find only indirect evidence for the role of nonthermal electrons in 
affecting the X-ray spectrum during stellar flares, this work bolsters conclusions from other 
authors that additional processes affect stellar X-ray spectra beyond the usual contributions 
of optically thin thermal emission in collisional ionization equilibrium. The X-ray to optical 
flux ratios seen in the early stages of the flare will hopefully be useful in constraining the 
contribution of flaring M stars to the population of unidentified transient hard X-ray sources. 

The frequency of occurrence of flares of this magnitude is unknown, but they may be an 
important factor in determining the habitability of planets around M dwarfs, as even a single 
occurrence of such a large radiative release can have disastrous consequences for terrestrial- 
like atmospheres on exoplanets orbiting M stars. As progress is made on identifying nearby 
stars hosting exoplanets, and because of the renewed emphasis on M dwarfs as host stars, 
such studies must be done in parallel with those attempting to understand the detailed 
nature and distribution as a function of luminosity of the stellar flares themselves. 

This work made use of data supplied by the UK Swift Science Data Centre at the Uni- 
versity of Leicester. The Konus-Wind experiment is supported by a Russian Space Agency 
contract and RFBR grant 09-02-00166a. We thank the Swift team for the ToO on Swift ~ 
one year after the flare to determine the quiescent level of EV Lac in several UVOT filters. 

Swift, Liverpool:2m EV Lac 
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Table 1. UV and Optical Magnitudes of EV Lac 



Filter Telescope 



Time Interval^ 



Value 



V Swift/UVOT 

white Swift/UVOT 



u Swift/UVOT 



uvwl Swift/UVOT 



uvm2 Swift/UVOT 



uvw2 Swift/UVOT 



LT 



T0+156.58:T0+164.78 7.2±0.2 

Q (UVOT Tool) 10.1 

T0+172.9:T0+174.3 < 7 

Q (UVOT Tool) 11.7 

T0+17670.5:T0+17970.3 12.96±0.01 

T0+17974.0:T0+18273.8 12.97±0.01 

T0+18277.6:T0+18577.4 12.95± 0.01 

T0+24289.8:T0+24589.6 12.89±0.01 

T0+24593.5:T0+24893.3 12.93±0.01 

T0+24897.0:T0+25145.9 12.93±0.01 

T0+36759.4:T0+37018.0 12.80±0.01 

Q (UVOT Tool) 13.2 

T0+13822.0:T0+14383.7 13.75±0.01 

T0+23383.6:T0+24283.3 13.59±0.01 

T0+35853.2:T0+36752.9 13.36±0.01 

Q (290 s)^ 14.03±0.01 

T0+12915.4:T0+13815.1 14.86±0.02 

T0+30983.6:T0+31047.9 14.76±0.05 

T0+34946.5:T0+35846.2 14.05±0.01 

T0+42546.2:T0+43008.0 14.77±0.02 

Q (334 s)^ 15.27±0.03 

T0+29164.3:T0+30064.1 14.90±0.01 

T0+40727.6:T0+41627.4 14.89±0.01 

Q (334 s)^ 15.10±0.02 

TO+900 (30s)^ 7.4±0.2 



'^Time in seconds since the BAT trigger 

'^Observations were made almost one year after the flare, to 
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determine the quiescent magnitudes in these filters. The time in 
parenthesis gives the exposure time for this observation. 

'^Observations at the Liverpool Telescope consisted of 3 10 sec- 
ond exposures taken roughly 15 minutes after the trigger. 
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Table 2. BAT Spectral Fit Results'' 



Time Interval^ 


Temp. 


VEM 


^2,c 


ix (14-: 


100 


keV) 


fx/fbol 


(s) 


(keV) 


(10^^ cm-3) 




(10-9 erg 


cm 


-2 S-1) 




TO-39: TO 


12.0it 


"•<J+-2.1 


46.87 


11 7+°-^ 






0.62 


TO: TO+112 


n.itll 


4 7+2-3 
^•'-1.4 


35.36 


^•-•-0.8 






0.45 


TO+112: TO+165.6 


8.9tli 


4.8+1^ 


24.45 


5 n+0-2 






0.26 


TO+165.6: TO+961 


10.0+^1 


0.81^^} 


34.91 


1 in+o-o'^ 






0.06 



^Global metallicity Z fixed to 0.4. 
'^Time in seconds since the BAT trigger 
■^37 degrees of freedom for each spectrum 



Table 3. XRT Spectral Fit Results 



Time Interval^ 
(s) 



kTi 
(keV) 



VEMi 



kT2 
(keV) 



VEM, 



(10 



54 pj^ 3^ 



X^d.o.f. 



fx (0.3-10 keV) 
erg cm~^ s~ 



(10-9 ... .^-2 .-!> 



T0+171:T0+306 

T0+306:T0+508 

T0+508:T0+671 

T0+671:T0+874 

T0+874:T0+1118 

T0+1118:T0+1451 

T0+1451:T0+1971 

T0+1971:T0+2215 

T0+2215:T0+2517 

T0+2517:T0+2798 

T0+6160:T0+4.3xl0^ 



h0.083 
-0.184 
h0.036 
-0.041 
h0.053 
-0.066 
h0.026 
-0.026 

U. (oy_Q025 

n 7c;q+o.024 

U. (O»_o.025 
hO.018 
-0.019 
[-0.021 
-0.022 
h0.021 
-0.022 
n 77Q+0.025 
"-"• ' ''-'-0.027 
n 797+0.025 
"-"• '^'-0.027 



1.196: 

0.990: 
0.956 
0.789 



0.772: 
0.759: 
0.771: 



0.58+° 

OMto 
0.32+° 
0.25+°: 
0.2471 
0.181° 

0-i6l°: 

0.0951 
0.0651 
0.0421 
0.0052 



,28 

18 

08 

,09 

,12 

08 

05 

04 

0.052 

0.008 

,04 

03 

03 

02 

0.02 

0.01 

0.010 

0.008 

0.007 

0.006 

+0.0009 

-0.0006 



Q OQO-^1.225 

0.000_Qg43 
7.06^ 



j+0.747 
'-0.634 
h0.781 
-0.656 
hO.445 
-0.303 

4.828l°;iti 

oqcri +0.221 
0.y01_Q 205 
hO.222 
-0.207 
[-0.147 
-0.139 
[-0.136 
-0.131 
h0.137 
-0.146 
h0.062 
-0.053 



6.036: 

5.222 



3.754: 
3.271: 
3.140: 

2.811: 

2.107 



-[-0.14 

-0.13 

+0.06 

-0.06 

0.07 

0.06 

+0.04 

-0.04 

+0.03 

-0.03 

+0.02 

-0.02 

+0.02 

-0.02 

+0.01 

-0.01 

+0.010 

-0.009 

2^8+°-°^° 



3.69 
2.35 
1.74 
1.51 
1.23 
0.91 
0.55 
0.40 
0.322 



477+0.110 
^■^' ' -0.099 

^62+°-^^^ 
n 442+0-129 

U-'±'±^-0.113 

h0.085 
-0.078 



0.439 



0.0096 



-0.0006 
-0.0006 



n s99+°-°'^6 

u.oyy„o.069 

U.^l'±„0.069 
+0.069 

-0.063 
+0.064 

-0.059 
+0.065 

-0.061 

h0.083 
-0.076 

h0.062 
-0.053 



0.429 
0.432 
0.497 
0.589 
0.382 



354.2,375 


24.51°:^ 


406.2,354 


15.81°:^ 


307.5,281 


11 1+0-2 

-'--'-•-'--0.3 


388.1,336 


9.2 ±0.2 


429.6,324 


7.4 ±0.1 


310.3,310 


5.2 ±0.1 


366.3,300 


3.3 ±0.1 


304.2,292 


2.22 ±0.03 


365.4,293 


1.77 ±0.03 


271.8,259 


1 Qq+0.02 
-'-•'-'^-0.03 


236.3,175 


n n^iT+o-ooi 

U.UO/^_0.001 



00 



^Time in seconds since the BAT trigger 
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Table 4. Time-Resolved Equivalent Widths of the Fe 6.4 keV line 



Time Interval^ 


Eb 


EW^ 


pb 


(«) 


(keV) 


(eV) 


(10~3 photons cm 


T0+171:T0+306 


6.44 ± 0.08 


146l?L 


272t\^l 


T0+306:T0+508 


6.421°:!} 


pro 9+64.8 
'J^-^-52.2 


6.96l|i 


T0+508:T0+671 


6.41 ±0.04 


320 ± 102 


26.41^:^ 


T0+671:T0+874 


(- ^9+0.07 


1031^^ 


6.5llt^? 


T0+874:T0+1118 


fi 4«+0-05 

0-46_o.o6 


1431}?^ 


6.86lto4 


T0+1118:T0+1451 


(- ^-1+0.15 


84.31^^3^ 


2.65111 


T0+1451:T0+1971 


"•^'-'-0.09 


Qll^f 


9 19+2.02 

^•-•-^-1.63 


T0+1971:T0+2215 


6.381°:}^ 


04 0+288.2 


1 20+°-^"^ 


T0+2215:T0+2517 


6.111°:?? 


001+52 
'-''-'-'--281.3 


1.65 ±0.85 


T0+2517:T0+2798 


6.22«:^^ 


2951^°^ 


0.951°:^! 



(%) (%) 



10.7/2 


99.672 


99.94 


1.8/2 


54.200 


55.42 


20.7/2 


99.994 


99.97 


5.2/2 


89.500 


99.19 


7.7/2 


94.563 


99.92 


3.5/2 


82.600 


97.69 


4.5/2 


84.100 


ND'^ 


3.4/2 


80.400 


ND'^ 


10.4/2 


98.502 


97.30 


4.8/2 


89.900 


ND'^ 



^Time in seconds since the BAT trigger 

^Errors are 3(7 

'^No Data - Monte Carlo calculations not done for this time interval due to low significance from 
F-test (see §3.2.2 for details) 
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Table 5. Spectral Fits to Time Interval T0+171:T0+961 s 



2T 



3T 



XRT+BAT 



kTi (keV) 






n Q8+0-02 


0.661°:°? 


VEMi (10^^ cm-3) 






n qc;+o-o5 


O-llTot 


kT2 (keV) 






6.68+°iI 


1 24+0-09 


VEMa (10^^ cm-3) 






2.10+^1 


Q OO+0.04 
"•'-'^-0.04 


kTs (kev) 






- 


7 Q4+0.31 


VEM3 (10^4 cm-3) 






- 


1 q«+o.05 

-■-•^0-0.04 


A 






r. 40+0.05 


n c;n+o-o5 


E (keV) 






6.43«:°^ 


6.43l°:i 


N (photons cm~^ s" 


-') 




(1.35l°;i)xl0-2 


(1.3ll°:i)xl0-^ 


X\ dof 






777.00, 589 


698.94, 587 


f(0.3-10 keV) (erg cm-^ g- 


-') 


(1.390+°-°??) X 10-8 


(1.395l°:r) xlO-8 


f( 14-1 00 keV) (erg 


cm ^ s 


-') 


(1.04«:°^)xlO-9 


(1.291°:°^) xio-9 


— XRT only — 


kTi (keV) 






n Qo+0.02 
U-9»_o.o3 


67+°-°^ 

U-O'-0.05 


VEMi (10^^ cm-3) 






n qc,+o-06 


12+0-06 


kT2 (keV) 






fi 67+0-59 
"•"'-0.28 


1 30+°-2 

-'-•'-"-'-0.08 


VEM2 (10^4 cm-3) 






9 1 n+0.03 
^•-'-'-'-0.03 


n Q7+0.05 
U-<J'-0.04 


kTg (kev) 






- 


8.221°:^^ 


VEM3 (10^4 cm-3) 






- 


1-921°:°^ 


A 






0.43±0.05 


0.50±0.05 


E (keV) 






6.431°:°! 


6.44±0.04 


N (photons cm"^ s" 


-') 




(1.35l°:i)xl0-^ 


(1.23l°:i)xl0-2 


X^ dof 






737.25, 550 


646.69, 548 


f(0.3-10 keV) (erg cm'^ s^ 


-') 


(1.3901°;°°^) xlO-« 


(1.410±0.008)xlO"S 
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Table 6. Flare Radiated Energy Estimates 



Energy Range 


Time Range^ 


Energy 


(keV) 


(s) 


(1033 erg) 


14-100 


T0-39:T0+961 


7.Q^ 


0.3-10 


T0-39:T0+165.5 


20^= 


0.3-10 


T0+171:T0+19792 


55^^ 


0.3-100 


T0+171:T0+961 


37^ 


0.3-10 


T0+171:T0+8300 


58f 



'^Time in seconds since the BAT trigger 

^Determined from fits to BAT spectra 

^Determined by extrapolating fits to BAT spec- 
tra in this energy range 

"^Determined from fits to XRT spectra 

'^Determined from joint fit to XRT and BAT 
spectra 

^Determined from exponential fit to flux vs 
time with a time constant of 647 s before 
T0-|-1000s, and a time constant of 1042 s from 
T0+1000s:T0+8300s 
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10^ 
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Fig. 1. — Light curves during the flare from EV Lac. Panels i-iii display emissions up to 3000 
s past the trigger time, TO of 2008-04-25T05:13:57, depicting from top to bottom, the flare 
response as measured by Konus/Wind in the 18-70 keV range (i), Swift/BAT from 14-30.1 
keV (a), Swift /XRT from 0.3-10 keV {in). Vertical dashed line indicates the BAT trigger 
time. The Konus/Wind data have been binned to 58.88 s intervals, and have been corrected 
for the time-of-flight between the Konus/Wind and Swift satellites: the horizontal dotted 
line indicates the background levels. Blue lines in the Swift/BAT and Swift/XRT panels 
indicate the times over which spectra were extracted. Panels iv and v display emissions 
up to 4x10^ s past TO. For the XRT light curve {iv), black points indicate data taken in 
Windowed Timing (WT) mode, red points indicate data taken in Photon Counting (PC) 
mode. The bottom panel {v) shows the UV and optical photometry for the flare. 
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Fig. 2. — Results of time-resolved spectroscopic analysis. Upper left panel details variation 
of plasma temperature as derived from BAT spectra (black circles) and the low and high 
temperatures determined from XRT spectra (red and blue circles, respectively). Horizontal 
bars give time range over which spectra were extracted; vertical bars are 90% error ranges. 
Upper right panel shows variation of flux in the 14-100 keV range (black circles) and 0.3-10 
keV range (red circles). Filled red circles correspond to measurements from XRT spectra; 
open red circles are extrapolations from the best fit model to the BAT data. Dashed line 
indicates bolometric flux of 1.7xl0~® erg cm~^ s~^. Lower left panel shows trend in volume 
emission measure derived from APEC model fits to BAT data (black circles), and the low 
and high temperature components fit to XRT spectra (red and blue circles, respectively). 
Lower right panel plots variation in metal abundance Z as determined from spectral fitting 
to XRT data. 
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Fig. 3. — Close-up of the Fe K region for the first five time intervals in the initial stages of 
the flare decay. Times are indicated in the top of the plot and refer to seconds since the 
BAT trigger time TO. Figure [H shows the XRT light curve with lines indicating the times 
when these spectra were extracted. Time intervals have been chosen so that spectra have a 
minimum of 20,000 counts per spectrum before background subtraction. Pluses are spectral 
data; dotted line is the best-fit two-temperature APEC model; dashed line is the Gaussian 
fit to the feature at 6.4 keV, and the solid histogram is the combination of the APEC and 
Gaussian models. 
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Fig. 4. — Modelling of 6.4 keV line flux for the first five time intervals during the flare 
decay; times are indicated corresponding to extracted spectra shown in Figure [31 Colored 
curves give expected run of Ka fluorescent line flux as a function of astrocentric angle for 
different flare loop heights. Dotted line indicates measured Ka flux and magenta region the 
3a uncertainty in the Ka flux. 
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Fig. 5. — Joint spectral fit to XRT and BAT data accumulated over the time interval 
T0+171:T0+961. A three temperature APEC model was used along with a Gaussian line 
at 6.4 keV. While a two temperature model was fit to the XRT spectra analyzed on finer 
time intervals, the addition of a third component was necessary. The highest temperature 
plasma indicated by fitting the XRT and BAT spectra simultaneously is the same as that 
from a fit to only the XRT spectrum, and there is no evidence for bremsstrahlung emission 
from a superhot plasma or nonthermal emission in the hard X-ray spectrum. 
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Fig. 6. — Variation of temperature versus square root of emission measure for times in the 
decay phase of the flare, using XRT spectra in the time range TO+171 s: TO+2798 s. Results 
have been derived from fitting XRT spectra to isothermal APEC models. Temporal evolution 
proceeds from the top right to the lower left. 
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Fig. 7. — A grid of parameter values for thick target bremsstrahlung emission was explored 
to determine if this was a viable mechanism to explain the excess Ka flux during the time 
interval TO+171 s:T0+961 s. Values of the total power in the electron beam Fq are plotted 
against the power-law index of the electron distribution function 6, for a cutoff energy Eq = 
100 keV, and column density in the Ka -emitting layer A^* =10^^. Plus symbols indicate 
values compatible with the excess Ka flux measured. Plus symbols over-plotted with open 
circles are grid points compatible with both the excess Ka flux and a 20 keV photon flux 
less than the amount 10^'^ photons cm^^ s^^ observed in the BAT spectrum. See §3.6 for 
more details. 
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Fig. 8. — Schematic showing the geometry between the flare photons, emitted as high en- 
ergy X-ray continuum, the region where the photospheric Ka photons are emitted, and the 
probable location of the white light flare. Different length scales used in the text are also 
indicated. The coronal loop is shown scaled to the stellar radius in approximate ratio to 
the values indicated from our analysis, and the foot-point area is a lower limit based on the 
white light flare emission. 



